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Of the various interesting solutions found in the 
two-dimensional complex Ginzburg-Landau equation for 
anisotropic systems, the phase-chaotic states show particu- 
larly novel features. They exist in a broader parameter range 
than in the isotropic case, and often even broader than in 
one dimension. They typically represent the global attractor 
of the system. There exist two variants of phase chaos: a 
quasi-one dimensional and a two-dimensional solution. The 
transition to defect chaos is of intermittent type. 

PACS: 47.54.-Hr, 05.45.-Fb, 47.20.Ky, 42.65.Sf 



The complex Ginzburg-Landau equation (CGLE) 
plays the role of a generalized normal form for spatially 
extended media in the vicinity of a supercritical Hopf bi- 
furcation involving a non-degenerate (oscillatory) mode. 
It has a wide range of applications extending from hy- 
drodynamic instabilities and nonlinear optics |^] 

to oscillatory chemical instabilities like the Belousov- 
Zhabotinsky reaction [Q or oxidation on catalytic sur- 
faces For a general review see e.g Q. 

The one-dimensional (ID) and the 2D isotropic cases 
have been investigated rather well |^]- |l6[. A number 
of results have also been obtained in 3D ]17| , p8[ |. Taking 
up some earlier work we recently reported about spi- 
rals and ordered defect chains in the anisotropic complex 
Ginzburg-Landau equation(ACGLE) H^] 



dtA =[! + (! + ibi)dl + (1 + ib2)dl - (1 + ic)\A\ 



]A. 
(1) 



Here A is the complex amplitude modulating the criti- 
cal mode in space and time. The usual reduced units are 
used. This equation was also studied in the context of de- 
fect chaos (DC) and wind-driven Eckmann boundary 
layers ||. 

Actually the range of applicability of Eq. (^ is consid- 
erable. The isotropic case, i.e. Eq.(^ with bi = b2 — b, 
can essentially be applied only to isotropic systems un- 
dergoing a spatially homogeneous Hopf bifurcation. A 
nonzero wavenumber Qc leads to traveling or standing 
waves, as in many hydrodynamic instabilities. Then, in 
systems that are isotropic in the basic state, one has a 
continuous degeneracy of the critical modes, which makes 



a more elaborate description necessary. In the presence 
of an anisotropy, like e.g. in the well-studied system of 
electro-convection in liquid crystals ||] , this degeneracy is 
typically lifted and Eq.(|l|) is appropriate. Also, of course, 
anisotropic systems with a Qc = bifurcation, as occur 
in oscillatory surface reactions require fei ^ 62- Tak- 
ing linear transformations of x and y into account the 
term involving second derivatives is general. Transform- 
ing into a co-moving frame a linear group velocity involv- 
ing a first space derivative vanishes. In the (common) 
situation of degeneracy between left- or right-traveling 
waves, we assume only one type to survive (which is of- 
ten the case). 

The ACGLE has a 2D wave-vector band of plane- 
wave solutions A — F X ex-pi{Qx + Py — ut), — 
1 - - P^.LU = c + {bi - c)Q^ + (62 - c)p2. They 
are stable against long-wavelength modulations when 
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)g2 + (1 + 2^1i^)p2 < I holds (generalized 



Eckhaus instability) while the Newell criterion 



l + b^c> 0, 
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is satisfied in both directions. From these relations one 
sees that the stable Q band shrinks to zero as 1 -I- 6ic — > 
0''" (Benjamin- Feir (BF) instability) with a similar behav- 
ior of the P band. Actually, the Eckhaus instability for 
{Q,P) ^ is of the convective type and plane waves can 
occur over a limited spatial extension in a larger range 

The bifurcation connected with this instability is su- 
percritical when one is at the BF limit or sufficiently 
near to it, i.e the amplitude of the destabilizing side- 
band modes actually saturates |2^. However, the re- 
sulting quasi-periodic solutions, as far as they are them- 
selves modulationally stable, have for vanishing (Q, P) a 
small basin of attraction in the BF unstable range, and in 
the studied ID and isotropic cases the relevant attractors 
turn out to be spatio-temporally chaotic. Nevertheless, 
since the amplitude \A\ saturates to a value near 1, only 
the phase $ (we write A = |^|expi<I>) is dynamically 
active. In ID the bifurcation at the BF instability, in- 
cluding slow modulations, is captured by the celebrated 
Kuramoto-Sivashinsky phase equation (see below). It ex- 
hibits the so-called phase chaos (PG) (or phase turbu- 
lence). 

PG in the ID CGLE was studied numerically first 
by Sakaguchi |Q, who also studied the breakdown and 
crossover to chaos involving phase slips (zeros of A in 
space-time) further away from the BF curve. This state 
is, in analogy to the 2D case (see below) , often referred to 
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as defect turbulence or defect chaos (DC) 0. The result- 
ing phase diagram was studied numerically in detail by 
Shraiman et al. ||2^, who discovered that for |6| > 1.8 the 
transition between PC and DC is continuous, whereas it 
is hysteretic with a bi-chaotic region otherwise (see Fig. 
|l|, dashed-dotted line and region marked bichaos ID). A 
detailed study with longer simulations and larger systems 
was performed by Egolf and Greenside [p6[. 
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FIG. 1. phase diagram for b2 = c 
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A rather exhaustive study in 2D (isotropic case bi = 
62) was presented by Manneville and Chate Here the 
region of PC is somewhat smaller than in ID (see Fig. 
|l|, broken line). Also, the transition is always hysteretic, 
which may be related to the fact that the zeros of A 
now correspond to topological defects. The breakdown 
of PC involves the creation of pairs of defects of opposite 
polarity which separate and loose correlation ("unbind"). 
Once initiated, the process is self sustaining leading to a 
nucleus and eventually to fronts that appear always to 
invade the PC state Thus in the isotropic case PC 
is never the globally stable attractor. 

Actually over much of the region where one has PC 
the global attractor is not DC as such, which appears 
only transiently, but rather a frozen state (vortex glass) 
with a disordered distribution of defects |^Jl^,|2^ . Every 
second defect emits a spiral wave of the type well known 
in the Eckhaus-stable range. The emitted waves remain 
intact over finite-sized cells by convective stabilization. 
Rotating spirals (time dependence oc exp icut) exist also 
in the ACGLE. In spirals the group velocity, which in 
plane waves is 2{bi — c)q in the x direction [2(62 — c)p 
in the y direction] is expected to point outward in all di- 
rections. In order to have coherent wavefronts one needs 
{hi — c)(&2 — c) > 0. Our simulations confirm that spirals 
are found only under this condition. Also the expected 
aspect ratio a/ (61 — c)/{b2 — c) of the equiphase lines of 
spirals is confirmed by the simulations. 

Our investigation was motivated in particular by the 
question of what happens in the parameter regime (bi — 



c)(&2 — c) < where spirals do not exist, and therefore 
also the existence of DC could be questionable. With 
this inequality the BF instability (necessary for PC) 
can only occur in one direction (we choose bic < — 1, 
i.e. instability in the x direction) and the anisotropy 
is "strong" |Q. Since the ACGLE has the symme- 
try (61,62,0,^) — > (— 61, — &2, — c, A*), we always chose 
&i > (in comparing with other works we transformed 
to this convention), and therefore &i — c > 0. 

The quick answer to the above question is actually 
quite simple: The system remains in PC "longer" than 
in the isotropic case, but eventually it does develop 
("anisotropic") DC. Since in DC defects actually hardly 
emit waves, in contrast to the situation in the vortex 
glass, no problem arises with opposite group velocities. 
The investigation led to surprises to be discussed now. 

We have performed detailed simulations of the ACGLE 
in systems of size L between 100 and 2500 dimension-less 
units with discretization At « 0.1 and Ax = L/N be- 
tween 0.3 and 5, where N is the number of Fourier modes 
in each direction of the pseudo-spectral algorithm used. 
We used periodic boundary conditions with initial condi- 
tions that imposed a zero phase difference across the sys- 
tem. Hence PC with a nonzero background wave-vector 
as studied recently in ID |14| was excluded. The results 
depend only weakly on the discretization and on system 
size (for sufficiently large systems). Choosing 61 ~ 62 
the results of ||2^ could be reproduced. Subsequently we 
changed 62 in the direction of c. This always increased 
the range of PC (i.e. |c| could be chosen larger). The 
limit of PC for the case 62 = c is depicted in Fig. |l|. To 
the left of the shaded region no defects were observed, 
to the right of it DC was found. The shaded region it- 
self is the parameter range where we found intermittency 
(see below). Note that for 61 > 2 even the ID limit of 
PC could be exceeded. Since it turns out that for the 
effect on PC the sign of 62 — c is, after all, not decisive, 
we in fact did many of the studies at 62 = c (< 0). A 
snapshot of the PC found there is shown in Fig. 0a. 




FIG. 2. Modulus \A\ for 61 = 5.0, 62 = c = -0.25, and 
system size L = 700 (white: !y4| = 1, black: \A\ « 0.9).a) : 
2D PC ("PCII") b): quasi-lD PC ("PCI"). 



We now come to the qualitative new features of 
anisotropic PC as extracted from our simulations per- 
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formed in the range 1 < 61 < 5 and 62 = c: 

(i) PC is the global attractor, i.e. with random initial 
conditions the system ends up in PC after the eventual 
annihilation of transient defects. This is in contrast to 
the isotropic case, where PC is never the global attractor. 

(ii) In the whole investigated range the transition be- 
tween PC and DC (as &2 = c is varied) goes through a 
stage of intermittency (Fig. ^ shaded region), which is 
not found in the isotropic case. In the intermittent state 
defect pairs are created in the form of bursts which subse- 
quently annihilate again, keeping the correlation between 
partners, i.e. defect pairs remain bounded. So in this 
regime, in spite of the presence of defects, phase coher- 
ence persists and the state should therefore be classified 
as PC. At a critical value of |c| (= c„) defects start to 
unbind rather fast, and this should be associated with 
the onset of DC. Recently a transition between two de- 
fect chaotic states in coupled Ginzburg-Landau equations 
was reported where one also sees this unbinding of pairs 
|§. _ 

In PC the spatial average of the amplitude |^| is very 
close to 1, see Fig. ||a (solid line). One finds a kink at 
the onset of intermittency from where on \A\ starts to 
drop faster. The limit of existence of PC can here be 
assigned to |c| slightly below 0.8. Also shown in Fig. ^ 
is the minimum |A|rrim of \A\ (broken line). Once |^| 
falls below |^| — 0.6 breakthrough to A = Q typically 
occurs. Figure ^ shows |v4|mi„ as a function of time in 
the intermittent range. 
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FIG. 3. a): \A\ (solid line) and lAjmin (broken line) de- 
crease with growing I&2I = |c| {L = 700, 61 = 4.0). b): time 
series of l^lmin for 61 — 2.0, b2 = c — —0.96 



(iii) In addition to the 2D PC discussed up to now 
("PCII") there exists close to the BF boundary and co- 
existent with PCII a strictly quasi-lD PC ("PCI") with 
spatial variations only in the unstable x direction (see 
Fig. ||b for a snapshot). It is obtained by initializing the 
system with a function A that differs from 1 only by small 
variations in x. PCI is only stable against small per- 
turbations in the y direction and easily transforms into 
PCII (it is metastable). At its limit of stability, which 
for bi = 5.0 is slightly below |c| = 0.3(= I62I), the trans- 
formation becomes spontaneous Because of CPU 
time limitations the transition could not be studied ex- 
tensively. 



Next we introduce a nonlinear phase equation which 
should yield a simplified description of phase chaos be- 
coming exact in the limit 1 -I- bic 0^. Using the stan- 
dard procedure Q one arrives at the following equation 
for the (strongly) anisotropic situation 

dt<s> = -\dM^ - D.4dt^ - gAd.^f + Dydl<f - (3) 

D^ = l + bic, Dy = l + b2C, D^4^bl{l + c^), 
9x^bi-c, gy^b2-c, a = bi{l + c^). 

The first three terms on the r.h.s. of Eq. ^ make 
up the ID Kuramoto-Sivashinski equation. The higher- 
order nonlinear terms proportional to a were included by 
Sakaguchi, who showed them to be responsible for the 
breakdown of PC, here implied by a blow up of the phase 
gradient Actually the last term in square brackets 
is formally of higher order than the others, but it could 
become important for small 61. In the stable y direction 
{Dy > 0!) it suffices to include the two terms shown, as 
done by Bar in the equation without the Sakaguchi terms 
Actually, in the parameter range studied by us, the 
term proportional to gy has little influence (for 62 = c it 
vanishes anyhow). 

By rescaling t, x, y, and <& one can scale the coefficients 
of the linear part and the term proportional to g^ to 1. 
Introducing the time scale r = D^i/lDxl'^ , the length 
scales become Ix = \/DxaJ\DA\ and ly = ^ DyDr^^|\Dx\^ 
which is supported by the simulations of PCII. Note that 
when the BF boundary is approached, where 0, ly 

diverges more rapidly than l^, so that PCII appears more 
and more one dimensional. Neglecting in Eq. (^ the 
last term in square brackets the only relevant parameter 
is the prefactor of the Sakaguchi terms, which becomes 

a = a\Dx\/{Doo4:gx)- 

Comparing PCII obtained from simulations of Eq. (|^) 
with that of the ACGLE we find satisfactory agreement 
except near to the breakdown (for not too large values of 
61). For bi — 2.0 and b2 — ewe find the breakdown of the 
phase description at 62 = —0.95 which is in fair agree- 
ment with the value found for the ACGLE. In a detailed 
study of the ID case at bi = 3.5 the authors found 
cid — —0.75 in the CGLE and —0.55 in the Sakaguchi 
equation, whereas we find C2d — —0.9 in the ACGLE 
(with &2 = c) compared to -0.75 with Eq. (|). PCI is 
also found in the phase equation and can at bi — 5.0 be 
maintained stably up to at least c = —0.28. The lowest- 
order description by Eq. (^ with only the first four terms 
on the r.h.s. has PCI and PCII as coexisting solutions. 

How can one understand the existence of PCI? For a 
stable ID solution, i.e. a solution with negative Lyapunov 
exponents, the (stable) existence of its quasi-lD analog is 
clear in the situation of a stable y direction (this is most 
easily seen in the phase equation). On the other hand. 
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in PCI one has positive Lyapunov exponents for fluc- 
tuations that vary only in the x direction, so there are 
also positive Lyapunov exponents for sufficiently small 
modulation wavenumber p in the y direction. However, 
this does not necessarily destroy PCI, since the only con- 
dition is that fluctuations flatten out in j/, even though 
they do not decay. We have confirmed by extensive simu- 
lations of PCI at 61 = 5.0, 62 = c = -0.26, L = 700, and 
N — 256 that small perturbations of the form exp ipy 
with Op < 0.2 (at p w 0.1) decay asymptotically in a dif- 
fusive manner with a phase diffusion constant around Dy. 
Under the same conditions stochastic perturbations (un- 
correlated on the discretized lattice in real space) decayed 
up to an amplitude < 0.01. Actually, one also expects 
solutions of Eq. (|l]) of the form A — exp {iPy)B{x,t) 
with phase-chaotic B to exist. Thus PCI presumably 
represents the center of a P band of phase- winding solu- 
tions. 

Finally we point out that the interpretation of the PCII 
<-> DC transition as a vortex binding-unbinding transi- 
tion probably allows to establish PCII as a thermody- 
namic phase that is qualitatively different from DC. In 
PCII, even if a defect pair is created, it remains bounded 
and annihilates again (the unbinding beyond c„ is a coop- 
erative phenomenon). The question of the conventional 
forms of PC representing such a state - in contrast to 
being just a (sometimes metastable) variant of DC with 
a very low rate of phase slips or defect pair creation - 
has indeed stimulated much of the previous research on 
PC Actually also PCI, although it appears to 

exist only metastably, can presumably be considered an 
independent thermodynamic phase because it differs in 
symmetry. 

Clearly much remains to be done. On one hand, find- 
ing criteria for the occurrence of PCI and methods to cal- 
culate the boundary of existence seems a most interesting 
problem. On the other hand, a detailed characterization 
of the PCII ^ DC transition appears desirable. 
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J. Neubauer, W. Pesch, and A. Rossberg. Extensive 
use of high-performance computer facilities at the LRZ 
Miinchen (Cray T90) and the HLRS Stuttgart (NEC 
SX4), as weU as financial support by DEC (Kr690/4) 
are gratefully acknowledged. 
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